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Abstract. By carrying out extensive lattice regularized diffusion Monte Carlo 
calculations, we study the spin and density dependence of the ground state energy 
for a quasi-one-dimensional electron gas, with harmonic transverse confinement 
and long-range 1/r interactions. We present a parametrization of the exchange- 
correlation energy suitable for spin density functional calculations, which fulfills 
exact low and high density limits. 
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i 1. Introduction 

. In this paper we present a parametrization for the exchange-correlation energy of a 

quasi-one-dimensional electron gas (IDEG) at arbitrary polarization. The electrons 
interact via a - potential and are confined to a line by a transverse harmonic potential 

, v(r±) = where b controls the thickness of the wire. Here and henceforth we use 

^ I the effective Bohr radius Oq = ^^^2 as unit of length and the effective Rydberg 

Ryd* = as unit of energy, where e is the dielectric constant of the embedding 
medium and to* is the effective electron mass. We can separate the transverse and 
' longitudinal parts of the Hamiltonian by assuming that the electrons are in the 

ground state of the 2D harmonic oscillator in the transverse direction. This is a 
good approximation provided that ^ where is the Wigner-Seitz radius. The 
above condition is met for low enough densities and thin enough wires. Thus, it is 
possible to integrate out the perpendicular degrees of freedom and work with a strictly 
one-dimensional Hamiltonian with electrons interacting via an effective potential given 

by 

F,(.) = #expf^')erfcfMV (1) 
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This potential has been widely used in previous works to model the IDEG, and we 
refer the reader to Ref. [T] for a detailed description of the Hamiltonian we study 
and the variational wave function we use. Here, we employed the lattice regularized 
diffusion Monte Carlo (LRDMC)[2] algorithm to compute the ground state energy of 
the system at different densities and polarizations, (C = {N^ — N-^)/N). In one 
dimension, this method provides the exact energy within the statistical accuracy, since 
the nodes of the ground state wave function are known exactly. 

Despite the huge amount of work done for ID systems with a 1/r interaction, 
[I1I31I11I51IS1I71IH] a spin density exchange-correlation functional is still lacking, and the 
Bethe ansatz solution is not available in this case. Fogler^ derived an approximate 
mapping of the problem with a realistic Coulomb interaction onto exactly solvable 
models of mathematical physics, but the relation is valid only for ultra-thin wires and 
requires a careful matching between different regimes. The quantum Monte Carlo 
framework can provide a parametrization valid in all regimes, but so far a functional 
has been derived only for an unpolarized wire.[T] The present work fills this gap, and 
we provide a spin dependent density functional for the exchange and correlation energy 
suitable for DFT calculations of these systems. Indeed, the DFT framework has been 
applied quite successfully in ID, [101 El 1121 HI] mainly on short-range ID problems 
where the homogeneous reference was known via Bethe ansatz. 

The paper is organized as follows. In Sec. [H we show the results for the ground 
state energy and give a parametrization for the exchange-correlation part, while in 
Sec. Owe present conclusions. In the Appendix we derive the polarization dependent 
random phase approximation (RPA) expression for the correlation energy, which is 
used to set the high density limit of our parametrization. 



2. Exchange-Correlation Energy and Construction of an LSDA Functional 

We study the ground state energy of the IDEG as a function of density and spin 
polarization, and find a parametrization for the exchange-correlation energy based 
on theoretically known properties of the electron gas in various limits. The best 
parameters for the exchange-correlation functional will be determined via a 
minimization of our LRDMC values for the total energy. 

Following the usual notation, we separate the total energy e into three parts: 

e(r„ C) = et{rsX) + ex(r., C) + ^drsX) (2) 

where Ct is the kinetic energy of the noninteracting system, e,j; is the exchange energy 
calculated for the noninteracting wavefunction and Cc is the correlation energy which 
includes corrections to both the potential energy and also the kinetic energy due to 
the interactions. The first two terms are known analytically, while the third one is 
fully determined by our numerical results for the total energy e. The kinetic energy 
reads 

7r^(l + 3C^) 
48rg ' 

while the exchange energy is 

F{x) = - / dyv{y)- 
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with v{x) — 2Ei{x'^) exp{x'^) the Fourier transform of the potential in ([T]), where Ei 
is the exponential integral function. 

To derive an accurate parametrization for the exchange-correlation energy, it is 
useful to study both the high and low density limits in order to include them in 
the actual functional. The high density limit is estimated with the random phase 
approximation (RPA), while the low density physics is obtained through a mapping 
onto an effective one dimensional Heisenberg model. 

The RPA is very successful in describing the energy of the homogeneous electron 
gas at high density [SI [T]. Here we present the main result valid for ^ 1 with the 
effective interaction in ([1]), while a detailed derivation is given in the [Appendix A] It 
is worth stressing that in the high density limit (small r^) the ID model with effective 
pair interactions given by (1) does not accurately describe electrons confined in a 
transverse harmonic potential, since the condition rs ^ ^ is in general manifestly 
violated, and the single subband approximation breaks down. The correlation energy 
evaluated within the RPA is 



where C = zi"^ {z)dz / {2'k%'^) w 4.9348/ (27r462). Though from JS]) the correlation 
energy may at first appear discontinuous at C = 1, ^^^^{I's, C) is in fact a continuous 
function of its variables as the two limiting behaviours in ([5]) clearly belong to different 
regions in the plane. 

The low density dependence of the correlation energy is difficult to determine 
since the effective coupling is very strong. This causes the electrons to repel each 
other and form a quasi- Wigner crystal. [3] As the exchange between the particles 
drops off very rapidly with the reduction in the density, different spin configurations 
become almost degenerate. However, the Lieb-Mattis theorem[T4 proves that in one 
dimension the ground state energy of a system of fermions corresponds to zero total 
spin. This theorem precludes the existence of a Bloch instability such as that predicted 
by an STLS-like theory. [5l [6] The low density spin dependence of the correlation 
energy can be determined approximately by noting that the spin sector of the IDEG 
can be mapped to that of a Heisenberg spin chain [15] with coupling J. In fact at 
these densities the electron gas is a quasi- Wigner crystal with local antiferromagnetic 
correlations. [16l lTll3| The Heisenberg coupling can be determined by an evaluation of 
the tunneling (exchange) rate between electrons via the WKB approximation, which 
gives an exponential suppression of J at low density as stated by the relation: (TTj 



where J* and v are interaction dependent constants. The energy dependence as 
function of J of the antiferromagnetic Heisenberg spin chain is known exactly from 
the Bethe ansatz.[T5] The difference in energies between the polarized and unpolarized 
spin chains turns out to be Jln2.[Tnj These relations define the spin dependence of 
the total energy of the electron gas at low density. Note that in order to provide the 
exponentially small spin dependence given by ((6|), the correlation energy must cancel 
the power law and logarithmic terms of both the exchange and kinetic terms. 

Our spin dependent exchange-correlation functional is built upon the 
parametrization of the exchange and correlation energy for the unpolarized (C — 0) 




(5) 
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and polarized (C = 1) wires, which reads: 



1 + d(;rs + e(^rj + /(^r^ 1 + h(;rf 

where the parameters are constrained to fulfill the high density limits of both exchange 
and correlation terms. Those limits imply the following conditions on the parameters: 

= , (8) 

,„=i±2±i!!lW2 (10) 

2 + 7 + 21n(7r6) 

^1= 2^2 (11) 

^ 2 + 7 + 2 ln(7rfo/2) , , , 

ci=-C/4+ ' ^.^ ^ U i+aiei, 13 

2 

50 = ^ with /3o > 1, (14) 

31 = ^ ■^it^ /3i > 1, (15) 

where 7 = 0.5772156649 is the Euler's constant. On the other hand, the large rs 
expansion of the expression in ([7]) goes as In r^/rs. Indeed, in previous work [SI [T] it 
was found that both the correlation and the exchange energies go as In rs/r^ at large 
Ts, with their ratio approaching a constant in that limit, a condition which is fulfilled 
by our parametrization. 

It is useful also to define a constrained exchange-correlation functional for the 
unpolarized case in a way that is accurate for low densities. Since the QMC 
calculations have lower variance for the fully polarized system, we define e^°"(rs,0) 
to be equal to the difference between the fully polarized and unpolarized energies of 
the antiferromagnetic spin chain with coupling constant J{rs), determined using the 
WKB approximation. Thus we rewrite the exchange-correlation functional for C = 
as 

4c"(-.,0) = ^^%^+(l ^) Le(r.,l)-J(r,)ln2-t--^), (16) 

where O and R are additional fitting parameters, and J{rs) is the same as in ^ with 
J* = 184.53 and ly — 2.84968 determined via the WKB approach for our potential in 
([T]). In this way both the high and the low density limits are fulfilled. 
Finally, the fully spin dependent density functional reads: 

e.e(r„ C) = elT{r,,0) + K{r,X) + c.(r., C) + ^ ^ ^t(r.)(i-\^\Y (1^) 
[2 ((1 - w{rM^ + wir^X") (e.c(r., 1) - el^^^ir^M - 2{K{r,,0 + c.(r„ C))] 
where the additional functions are: 



b^r, + CQrl 9crs In 
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w{rs) = 
Cz{r s,C) 




e 



(19) 
(20) 



(21) 



(22) 



Cz{fs,C) is the small C expansion of the correlation energy around C = 0, while 
hz {^8 , C) is the variation of the exchange energy with respect to ( = at fixed 
Ts- Both expressions are taken in the high density limit. hz{rs,C) includes another 
parametric function (/icorr(''"s, C)) which accounts for the non analytic behaviour of 
the exchange energy around C = 1 at = 0. The form of (I17p was chosen to constrain 
the parametrization to attain energies determined by (fT6|) and This allows the 
parametrization to in principle satisfy the non analytic behaviour of the correlation 
energy at high density and C = li while the low density behaviour is fulfilled by the 
mapping onto the Heisenberg model. Even if the parametrization looks complex at 
the first glance, there are only 21 independent parameters. 

We have carried out extensive LRDMC simulations to find the best fitting 
parameters for our parametrization. We note that there is another "external" 
parameter 6, which sets the effective thickness of the wire and therefore defines 
the interparticle potential. It is of course possible to derive the parametrization for 
different widths, but here we chose to work with 6=1, which is close to the usual 
thickness of wires realized in semiconductor nanodevices.[20] The calculations for = 1 
yield a series of total energies as a function of density and spin polarization. 

Great care is taken to remove all biases in the LRDMC calculations of the energy. 
The lattice space error is removed by calculating the energy for different lattice 
spacings and extrapolating the results with a quadratic fit in the lattice space. Finite 
size effects are removed by calculating the energy for several numbers of electrons and 
extrapolating the result to the thermodynamic limit by fitting the data to the form 



where E is the energy extrapolated to the thermodynamic limit, N is the number 
of electrons in the calculation, and the constants ci and C2 are fitting parameters 
determining the size of the one-body and two-body finite size corrections. Additionally 
the number of electrons N is chosen in each calculation so that the number of electrons 
in each spin species is odd, thus avoiding degeneracy effects. The form in ([23|) is 
obtained by following the finite size analysis described in Ref. |21j . 

Our results are plotted in Fig. [1] and Fig. [2] which show the behaviour of the 
correlation energy as a function of the density and the polarization respectively. The 
correlation energy at high density (r^ = 0.1) as a function of the polarization shows 
vestiges of the non-analyticity in the correlation energy at ^ = 1 for ^ (see ([5])). 

Tables [1] and [2] present the various parameters that are obtained by a least- 
squares minimization fitting of the LRDMC values for the exchange-correlation energy 
computed at 17 different densities ranging from — 0.1 to = 50. From rg = 0.1 to 
Ts = 1.5 nine values of the polarization were used equally spaced from ^ = to ( = 1. 
For rg > 1.5, five polarizations C = 0, ;|,^,|, and 1 were used. These parameters 
produce a fit that has a reduced of 5.3 and an overall accuracy on the order of 
10~^i?yd*. The exchange correlation energy is plotted at several values of the density 
in Fig. [3] As one can see, it is in a good agreement with the parametrization at all 
densities. 



E{N) + 



C2\/lnN ci 



(23) 
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Figure 1. The correlation energy of the electron gas as a function of the density 
rs is plotted for b = 1 at five values of the polarization, ^. 
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Figure 2. The correlation energy of the electron g function of the 

polarization ( is plotted for 6 = 1 at four values of the density, r^. The upper 
right panel shows evidence of the RPA behaviour, quadratically as a function of 
f near C = and then assuming an abrupt change around C = 1- The upper right 
and lower left panels plot the intermediate case, while the lower right panel shows 
the f dependence at low density. 



3. Conclusions 

In this paper we have presented results for properties of a quasi-one-dimensional 
electron gas, with harmonic transverse confinement and long-range 1/r interactions, 
which is a model for confined semiconductor structures. By carrying out extensive 
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Table 1. Parameter Values for the Fit of e!^°"(rs,0) and excirs, 1) 
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Table 2. Other parameters of the parametrization 





Figure 3. Exchange-correlation energy exc vs the polarization ( at various 
densities. The solid line comes from the parametrization while the points come 
from QMC calculations. Their error bars are smaller than the point size. 



lattice regularized diffusion Monte Carlo calculations, we have determined the ground 
state energy as a function of spin and density, and we have presented a parameterized 
fit to the Monte Carlo data that can be used as a local density functional for exchange 
and correlation in spin density functional calculations. The form is given in ()17II22|) . 
It fulfills the high density limits of both exchange and correlation energies around 
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^ = and at C = 1. At low density a mapping to an Heisenberg spin chain has been 
used to work out the C dependence, while the determination of J comes from WKB 
calculations. The parameters are given in Tables [1] and [2] for a typical wire width 
6=1. The overall accuracy of the fit is on the order of lO^^Ryd*. 
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Appendix A. RPA calculation of the spin dependent correlation energy 

In this appendix we compute the correlation energy of a spin polarized IDEG in the 
high density limit, using the random phase approximation (RPA). We start from the 
general expression of the RPA correlation energy |22j: 



^RPA 



= — / dk 

J-oc 



1 la f+°° 

e{k)^ —— dX ln{l-i{kb)x°{k,ikX)) + ij{kb)x°{k,ik\), (A.l) 

where v{kb) is the Fourier transform of the potential, and X° = + Xj* is the real 
part of the density-density response function for the free ID electron gas: 

with Vp the Fermi velocity of the a (=t,i) component. After some algebra, and a 
change of variables {k = kpq,u) = ikpqvpu), (jA.l|) can be rewritten at the leading 
order as follows: 

with a = A^j-K in ID. The derivation follows the work of Gell-Mann and Brueckner [53] 
in 3D, and Rajagopal and Kimball [24 in 2D. The "propagator" (5^(m) depends now 
on the spin polarization, and reads: 

Q^(m)=/ dk \ dt Uk){l- U{k + q))e~'*^''exp{-\t\{-q^ + kq)), (A.4) 



— oo 



where /|(x) = 9{\x\ — (1 + C)): fl{x) = d{\x\ — (1 — C)) are the zero temperature Fermi 
distributions for the two spin components, 9 being the step function: 

(a.5) 

In order to factor out explicitly the order dependence in (|A.3p . we apply another 
change of variables, by rescaling q {q ^ ^F^?)' ^^'^ integrate over u. After these 
steps, the RPA correlation energy reads: 
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where we have defined the set of functions: 

Q J —oc J —oo 

q^ + qik,+k2)' ^^-^^ 
where the C dependence is included in the zero temperature Fermi distributions fa{k)- 
From the above equation it is apparent that = F^.[- 

For C = 1, F^^^{q,l) ^ 0, while F^^^[q,l) = F^^^{q^) = Vq. Since 
-F|.|(0, 1) = 7r/2, for the fully polarized IDEG we obtain: 

efP\r^X=l) = -^rl (A.8) 

a result which is in agreement with the mean spherical approximation. [S] 

To evaluate E^^^{rs, C) a-t intermediate polarizations, we need to compute the 
limits: 

lim Fa.a' {x, with C < 1. (A.9) 

It turns out that Ft,t(0,C) = 7r/(l + (), -FT,i(0,C) - ^, and ^,^(0,0 - V(l - 0- 
Thus, our final result for the spin dependent RPA correlation energy is the following: 

A n I 1 ^„2 



er%.,0 = \ + ifr <<2{1^0b/a ^^^^^ 

I Stri-b^ if C - 1 

Notice that when ^ = we recover the RPA correlation energy for the unpolarized 
IDEG derived in Ref. [1]. 
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